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ABSTRACT 


A  method  to  reduce  the  model  error  of  the  Hill  equation  for 
long  relative  distances  is  proposed  in  this  paper;  this  is  done 
by  coordinate  shifting  to  reduce  the  long  relative  distance  of 
the  original  system  to  a  small  one  in  a  new  system.  That  is,  as 
the  relative  distance  of  the  interceptor  reaches  a  particular 
precision  parameter  value  pj  (pj  is  defined  as  a  controlled 
precision  parameter,  for  example,  Pj=500km) ,  an  imaginary  moving 
reference  system  with  its  origin  at  the  interceptor's  initial 
position  is  set  up,  and  the  relative  motion  parameters  of  the 
interceptor  in  its  original  system  are  converted  into  those  in 
the  new  system  (obviously  the  initial  relative  distance  in  the 
new  system  will  be  .  As  a  result,  the  long  distance 


1 


problem  becomes  a  small  one,  and  the  descriptive  precision  of  the 
Hill  equation  is  improved.  This  paper  provides  a  theoretical 
analysis  for  this  method  and  two  examples  of  computer  simulation. 
This  method  does  indeed  reduce  model  error  for  long  relative 
distances  in  the  Hill  equation,  producing  satisfactory  results. 

Key  words:  space  interception,  satellite  orbit,  relative  motion. 
Hill  equation.  Type  Number  TJ861 

With  the  development  of  modern  aeronautic  technology, 
spacecraft  must  engage  in  a  variety  of  motion  and  flight  duties, 
such  as  fixed  time  target  interception.  The  problem  of  fixed 
time  interception  can  be  studied  in  inertial  reference  systems 
([1]  introduces  the  classic  Herrick  method,  Godal  method,  and 
Lambert  flight  time  theorem,  etc.).  This  problem  can  also  be 
studied  from  an  appropriate  motion  reference  system,  and 
reference  [1]  offered  a  description  of  two  Hill  equations  for 
relative  spacecraft  motion. 

Inertial  motion  system  research  is  based  on  orbital  dynamics 
theory,  and  a  precise  solution  for  the  fixed  intercept  time  Cit  is 
found  through  an  integrated  motion  equation  or  an  iterative 
solution  of  the  algebraic  transcendental  equation.  One  cannot 
directly  write  out  the  analytical  solution's  results,  so  it  is 
inconvenient  to  analyze  the  problem.  Motion  systems  in  the  Hill 
equation  make  it  possible  to  derive  a  relative  motion  equation 
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through  gravitational  linearization;  it  is  superior  in  that  it 
is  easy  to  derive  a  solution  which  closely  approximates  the 
analytical  solution.  But  it  only  has  a  highly  precise  solution 
for  short  relative  distances.  When  the  relative  distance  is 
longer  {for  example,  when  p>1000km) ,  the  model  error  in  the  Hill 
equation  is  relatively  large,  and  the  precision  of  the  solution 
is  decreased. 

In  actual  engineering  problems,  sooner  or  later  we  always 
encounter  situations  where  spacecraft  flight  time  {and,  therefore 
relative  distance)  is  relatively  long,  such  as  in  the  case  of 
long-distance  interception.  In  the  initial  study  of  the  problem, 
we  hope  to  have  a  fast  and  convenient  methodology  such  as  the 
Hill  equation,  but  if  we  then  use  the  simple  Hill  equation,  the 
error  results  caused  by  model  error  are  intolerable. 

The  critical  point  in  reducing  model  error  is  in  whether  or 
not  we  can  make  a  long  relative  distance  problem  into  a  small 
relative  distance  problem  through  the  use  of  some  type  of  linear 
conversion  {single  or  multiple) .  Assuming  the  interceptor  is 
released  from  a  satellite,  after  time  tj  its  relative  distance 
P=Pl  (Pi  being  a  precise  control  parameter,  for  example 
P2=500km) ,  and  constructing  an  imaginary  satellite  target 
trajectory  coordinate  system  {motion  reference  system)  at  this 
point  and  at  this  time,  then  at  this  specific  moment  the 
interceptor  is  at  relative  position  P()=0  in  this  new  system. 
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Size  of  the  interceptor's  distance,  we  can  always  make  the  long 
relative  distance  into  a  limited  small  relative  distance  through 
a  limited  frequency  of  this  kind  of  target  coordinate  system 
conversion,  so  that  every  time  we  use  the  Hill  equation  we  can 
guarantee  p<pp  and  reduce  its  model  error. 

1.  The  Hill  Equation 


In  a  distance-squared  inverted  gravitational  field,  the 
satellite  revolves  around  trajectory  I.  The  satellite  orbit 
system  0-XYZ  is  the  target  motion  coordinate  system.  If  the 
spacecraft  is  moving  close  to  the  satellite,  then  next  to  its 
orbital  path  there  is  a  small  bias  error  movement.  At  time  t  the 
geocentric  vector  is  r2,  and  the  satellite  location  vector  rj  is 
relative  to  the  satellite  location  vector  ^=r2-ry  This  is  shown 
in  Fig.  1.  If  the  spacecraft  and  the  satellite  are  both  subject 
to  the  function  of  the  distance  square  inverted  central  force 
gravitational  field,  then: 


dt*  + 


If  (i)  and  0  are  the  motion  coordinate  rotation  angle  of 
acceleration  and  angle  of  speed,  we  can  obtain  the  vector  form 
relative  motion  equation  as: 


^  +  2 
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with  target  coordinates  {X,  Y,  Z)^  in  the  spacecraft  motion 
system  0-XYZ,  we  can  use  gravitational  error  linearization  to 
obtain  the  spacecraft  relative  motion  differential  equation 
(assuming  p<<rj)  : 

jX  =  3n*X  +  2«Y 

jv=-2«X  (1) 

l2=-«*Z 


n  in  the  equation  is  the  trajectory  path  angular  speed  of  the 
satellite,  986005 XI 0‘«  m’/s* 

is  the  gravitational  constant,  and  M  is  earth's  mass.  This  group 
of  constant  coefficient  differential  equations  for  linear 
integration  approximates  a  description  of  the  relative  motion  of 
two  spacecraft,  and  is  called  the  Hill  equation.  The  smaller  the 
relative  distance  p,  the  more  precise  the  description. 

making  ?  =  (X.Y.Z.)*',  V  =  (X.Y.Z)’' 


Integrating  Eq.  (1)  yields  the  spacecraft  relative  motion 
parameter  at  time  t: 


Pr  ' 

Pw 

Vr. 

^-(r)J 

r<»>. 

pjQj  and  V|Qj  in  the  equation  are  the  initial  parameters  for  the 
spacecraft  at  time  t=0,  and  are  the  relative  motion 
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condition  shift  matrices. 


Kir)  = 


4  —  ZcosnT  0  0 

6(sin7ir  ~  nr')  1  0 

0  0  cosnt 

2(1  “  cosTzr) 


smnr 

n 

2(1  —  zosnz) 


n 

Asinnv 

n 


-  3r 


sinnr 


K^r)  = 


Kir)  = 


Snsinnr  0  0 

6n(cos«r  —  1)  0  0 

0  O'  —  nsinnr 

cosnr  2sinnr  0 

—  2sinnr  4cosnr  —  3  0 

0  0  cosnr 


Through  the  linearization  of  gravitational  error,  the  Hill 
equation  yields  an  analytical  solution.  From  Eq.  (1)  we  can  see 
that  relative  motion  possesses  the  following  special 
characteristics:  from  it  we  can  separate  solutions  for  two 
mutually  independent  movements,  the  longitudinal  plane  (XOY)  and 
lateral  (for  direction  Z) . 

2.  Target  Coordinate  Conversion  and  Relative  Parameter  Conversion 
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To  simplify  the  problem,  this  paper  will  only  study  the 
longitudinal  problem  at  present  and  not  consider  the  lateral 
(direction  Z)  movement,  and  continue  to  present  an  analysis  based 
on  a  three-dimensional  vector  form  (Z  direction  value  is  assumed 
at  0)  . 


Suppose  the  satellite  is  in  motion  in  orbit  1,  and  we 
release  a  flight  craft  with  relative  speed  of  The  motion 

system  at  time  t=0  is  0-XYZ,  and  the  satellite  latitude  is  Ug. 
After  time  Tj,  0-XYZ  follows  the  satellite  to  fly  to  position 
O'-X'Y'Z'.  At  this  time,  the  traveling  craft's  relative  position 
vector  is  and  let  us  assume  that  at  this  time  p=pj,  as  shown 
in  Fig.  2. 


Fig.  1.  Satellite  trajectory  Fig.  2.  Target  coordinate 
system  conversion 


At  Oj,  we  construct  an  imaginary  motion  reference  system 
Oj-XjY^Zj ,  and  we  can  easily  know 
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rx  =  Re H 


/^xO,X'  =  nsjTi 

^  —  2r,/)cosfl 

Re  in  the  equation  is  the  earth's  radius,  and  H  is  the  orbital 
height  of  the  satellite. 

Mark  that  at  time  Tj,  the  traveling  craft's  relative  motion 
parameters  are  p,  and  V  in  O'-X'Y'Z',  then  we  have 

?(rl)  = 


Assume  that  the  three  components  of  in  O'-X'Y'Z'  are 

«  =  y  +  arctgI>(M)//>(M)] 


du  =  sign[/0(j.,)]arccos 


A  +  A- 


2r,r, 


duu  =  d«  +  ns\^i 


Below,  we  convert  the  traveling  craft's  relative  parameters 
^  .  .  . 

p  and  V  in  the  original  system  to  the  relative  parameters 

PliO)  and  VjiDj  in  the  new  system.  Obviously,  Pj(QfO  and  Vj(u)  can  be 

derived  as  follows: 
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Take  Oe-  5tiC  as  an  inertial  system,  from 


=  frl  +  =^rt+  ^.2 

to  obtain 


ns\ri 

0 

0 

”52^2 

0 


+  -B(-0  +  «5,rj)^rl 


1(0) 


(3) 


In  the  equation 


nsz 


-  "J f^/r\ 


B 


idm'^ 


cosdu 

sindtf 

0 


—  sindw  0 
cosdu  0 
0  1 


®(du)  is  the 


Oi-XiYiZi 


target  coordinate  shift  matrix  between  O'-X'Y'Z' 


and 


From  Eq.  (3) ,  we  simplify  to  obtain 


^UO)  “ 


0 

(njj  — 
0 


(4) 


This  is  the  method  of  expressing  the  relative  speed  in  the  new 
system. 


3.  Two  Functional  Examples 
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We  use  a  satellite  released  interceptor  terminating  a 
missile  in  an  unpowered  missile  trajectory  as  an  example.  The 
satellite  orbit  path  is  a  polar  orbit,  at  a  height  of  H=640km. 
Time  of  interceptor  release  is  t=0,  latitude  is  Uq=50®. 

The  missile  is  launched  due  north,  and  its  termination  point 
motion  parameters  are  Hj=170km,  speed  V^=7390m/s,  local-speed 
angle  of  inclination  0^20.6®,  and  latitude  0j^=9O°.  Set  the 
interceptor's  largest  attack  impact  value  at  6421m/s,  and  the 
longest  interception  time  at  1500s.  The  satellite  orbit  path 
plane  and  missile  launch  plane  match,  and  perform  a 
unidirectional  chase  intercept,  without  calculating  earth 
flatness,  or  considering  earth  rotation.  This  paper  offers  three 
sets  of  computer  simulation  results:  the  first  set  is  based  on  a 
direct  calculation  of  the  Hill  equation  (indicated  as  H) ,  the 
second  uses  the  calculated  results  from  an  inertial  system  as  a 
standard  parameter  (indicated  as  HO) ,  and  the  third  offers  the 
method  in  this  paper  for  calculation  results  (indicated  as  HA) , 
and  graphs  the  deviation  curve  of  HO  with  respect  to  H  and  HA. 


Example  1  sets  intercept  time  at  t=1500s  and  uses  the  above- 
mentioned  three  types  of  methods  to  find  the  intercept  target's 
needed  attack  impact  speed  increase  Av.  Of  these,  the  HA 
simulation  process  precision  control  parameter  p2=200km,  and  the 
target  coordinate  system  shift  frequency  is  n  (times).  The 
results  are  as  follows: 
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HO.AV  =  (-  1506.77,1022.  12,0)’^ 


H,  AV  =  (-  1448.03,1062.  59,0)^ 

Note  that  5V^  is  the  H  deviation,  57^=71 . 33in/s ,  SVj  is  the  HA 
deviation,  the  5Vj,  5V2~P2  and  n^j  curves  are  as  shown  in  Figs.  3 
and  4. 


Example  2  sets  the  intercept  attack  impact  speed  increase  as 
follows:  Av=  (-1000 , 1000 , 0 ) '*'  m/s,  obtaining  pj=200km.  The 
intercept  termination  point  target  coordinate  (error  distance) 
follows  the  change  in  time  t  (and  relative  distance  p)  as  shown 
Fig.  5  (H,  HA)  .  The  HA  coordinate  shift  freguency  curve  n~p 
is  shown  in  Fig.  6. 


Fig.  3.  curve  Fig.  4.  n/sp^  curve 
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p(xiO’km)  p(xiO’km) 

Fig.  5.  Ap^  curve  Fig.  6.  ri/^  curve 


4.  Conclusions 

(1)  HA  can  effectively  reduce  the  model  error  in  the  Hill 
equation  for  long  relative  distances,  and  increase  its 
descriptive  precision.  Figs.  3  and  5  show  satisfactory  results. 

(2)  The  precision  control  parameters  for  HA  can  generally 
be  100«^00km,  and  at  this  time  can  guarantee  adequate  precision 
of  the  Hill  equation  for  long-distance  problems. 

(3)  HA  does  not  need  numerical  differentials  or  iterated 
solution  algebraic  transcendental  equations;  one  only  needs  to 
borrow  the  Hill  equation  for  an  analytical  solution  and  limited 
frequency  linear  conversion,  and  does  not  need  to  consider  the 
type  of  interceptor  trajectory  curve,  making  it  easy  to  program 
in-flight  processing. 

(4)  HA  has  a  definite  practical  application  value,  and  can 
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be  broadly  used  for  problems  such  as  interception  targets , 
non-orbital  motion,  orbital  shift,  and  orbital  path  intersection, 
etc . 


I  hereby  extend  my  deepest  gratitude  to  Professor  Dong  for 
his  assistance  in  the  research  for  this  paper. 
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